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I .  INTRODUCTION 

Computer-to-computer  communications  characterizes  the 
fundamental  link  for  message  and  data  traffic  in  many  modern 
military  communications  scenarios.  In  some  cases  the  link  is 
one  of  many  in  an  extensive  network,  in  others  it  may 
represent  a  single  point-to-point  communications  link. 
Regardless,  the  objective  is  to  transmit  bits  with  a 
satisfactorily  low  rate  of  errors  from  one  terminal  to 
another. 

The  medium  for  transmitting  the  data  may  be  wire  or 
optical  fiber,  that  is  a  low  pass  channel,  or  it  may  be 
radio  frequency  or  acoustic,  that  is  a  bandpass  channel.  In 
either  case,  the  actual  signal  used  to  carry  the  data  bits 
must  be  a  properly  modulated  analog  signal  with  sufficient 
energy  and  positioned  in  the  frequency  spectrum  to  propagate 
effectively  through  the  available  medium.  In  RF  links  the 
signal  will  be  further  restricted  to  an  assigned  frequency 
slot  or  channel  in  the  spectrum.  For  network  links 
transmission  may  be  restricted  to  assigned  time  slots. 

What  is  needed  is  a  signal  modulation  format  that  is 
readily  adaptable  to  a  variety  of  link  environments,  does 
not  require  special  purpose  MODEMS  to  translate  between  the 
digital  and  analog  domains,  can  emulate  most  existing  signal 
modulation  formats  as  well  as  generate  entirely  new  formats 


and  whose  descriptive  language  is  the  language  of  Digital 
Signal  Processing  (DSP) .  Multi-Frequency  Modulation  (MFM)  is 
such  a  technique  and   research  into  its  implementation  as  a 
practical  communications  method  is  the  subject  of  this 
report. 

Classical  modulation  methods  for  bandpass  channels  use 
amplitude  and/or  phase  to  carry  signal  information  on  a 
carrier  wave  in  the  channel  band.  When  the  information 
source  is  a  finite  alphabet,  as  with  data  or  quantized 
analog  sources  such  as  digitized  speech  or  video,  then  only 
a  finite  number  of  signal  states  are  needed  to  represent,  or 
code  the  source. (The  presence  of  one  or  the  other  of  two 
frequencies,  FSK,  is  also  used  to  transmit  binary  data,  but 
this  is  less  efficient  than  phase  shift  keying,  PSK, [ref 

1].) 

In  order  to  accommodate  multiple  sources,  individual 
channels  are  multiplexed  together  either  in  the  frequency 
domain, (FDM)  or  time  domain, (TDM)  or  both.  With  multiple 
users  on  a  network,  Time  Domain  Multiple  Access,  (TDMA) ,  or 
Frequency  Domain  Multiple  Access,  (FDMA)  is  used  to  allow 
each  data  source  the  opportunity  to  transmit  on  the  net. 

Modulation  and  multiplexing  equipment  for  communications 
links  carrying  many  channels  is  primarily  analog  circuitry 
and  is  very  expensive  to  purchase  and  maintain.  Both 
commercial  and  military  communications  system  developers 


desire  to  make  their  equipment  as  digitally  oriented  as 
possible  in  order  to  increase  its  flexibility  and  to 
automate  its  testing  [ref  2].  This  trend  is  being  followed 
at  SATCOM  earth  stations,  at  military  switching  centers, 
NAVCOMSTA ' s  for  example,  on  ships  and  aboard  aircraft. 

It  is  our  view  that  ultimate  digitization  is  achieved  by 
having  the  host  computers  absorb  the  functions  of  modulation 
and  multiplexing,  demultiplexing  and  demodulation.  MFM  is  a 
modulation  technique  that  is  ideally  suited  to  this  task 
because  its  basic  structure  is  one  of  time  and  frequency 
slots.  In  MFM  the  signals  are  directly  encoded  and 
modulated,  and  decoded  and  demodulated  using  digital  signal 
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Figure  1  MFM  Signal  Packet 


processing  techniques  within  the   host  computers. 
Connections  to  the  analog  link  are  made  thru  D/A  and  A/D 
chips  on  plug-in  boards  over  DMA  I/O  channels.   Over  eighty 
percent  of  current  digital  communications  signals,  including 
spread-spectrum  signals,  time  and  frequency  slotted  signals, 
plus  many  entirely  new  signal  formats,  are  simply  generated 
and  demodulated  using  MFM. 

In  MFM,  signal  "packets"  arbitrarily  located  in  the 
frequency  spectrum  and  in  time  are  created  as  shown  in 
Fig.l.  Each  packet  consists  of  L  bauds  of  K  tones.  These  LK 
subsignals  form  an  orthogonal  signal  set.  Each  subsignal  may 
be  independently  modulated  with  phase  and  amplitude 
information. 

From  a  theoretical  standpoint,  recent  research  has  shown 
that  multitone  signals  are  more  efficient  than  single 
carrier  signals  of  the  same  total  bandwidth  in  a  general 
linear  channel,  particuarily  if  there  is  a  deep  null  in  the 
frequency  response  [ref  3].  From  a  practical  point  of  view, 
we  have  suggested  that  MFQPSK  signals  are  ideally  suited 
for  use  in  high  data  rate  acoustic  burst  communications  from 
moving  platforms  using  medium  frequency,  relatively  wideband 
(25%  bandwidth)  links  [ref  4].  The  appealing  thing  about  MFM 
signals  is  the  simple  way  they  can  be  generated, 
demodulated,  and  equalized  using  standard  digital  signal 
processing  algorithims  which  are  now  available  for  nearly 
all  industry  standard  digital  computers  in  a  variety  of  high 


level  languages.  Real  time  processing  at  high  continuous 
data  rates  is  still  difficult,  however  improvements  in  this 
area  are  inevitable.  A  variety  of  array  processor  cards  are 
available  now  as  add-on  cards  and  new  designs  (the  NEXT 
computer  for  example)  are  incorporating  DSP  processors  into 
the  computer's  basic  architecture.  Thus,  although  current 
technology  limits  direct  bandpass  generation  of  MFM  signals 
to  the  tens  to  perhaps  several  hundereds  of  kilohertz  range, 
generation  of  wideband  HF  radio  signals  (  signals  in  the 
tens  of  MHz  range)  should  be  quite  feasible  within  the 
decade.  Of  course,  the  audio  range  signals  we  can  generate 
now  are  useful  in  acoustic  communications,  an  application 
that  is  being  pursued  by  us  in  conjuction  with  the  Naval 
Ocean  Systems  Center  in  San  Diego. 

II.  THEORY  OF  MULTI -FREQUENCY  MODULATION 

Generation  and  Demodulation  of  MFM 

The  following  definitions  are  used  in  MFM  (Please  refer 
to  Figure  1) : 

T:   Packet  length  in  seconds 

AT:  Baud  length  in  seconds 

kx:  Baud  length  in  number  of  samples 

L:  Number  of  bauds  per  packet 

At:  Time  between  samples  in  seconds 


fx=l/At:  Sampling  or  clock  frequency  for  D/A  and  A/D 
conversion  in  Hz. 

Af=l/AT:  Frequency  spacing  (minimum)  between  MFM  tones. 
K:  Number  of  MFM  tones. 

Since  At=AT/kx,  the  sample  frequency  f^J^Af. 
Consequently,  there  are  a  maximum  of  kx/2  tones  spaced  at 
intervals  of  Af  Hz  between  dc  and  less  than  fJ2,    the 
Nyquist  frequency,  that  can  carry  amplitude  and  phase 
information  during  each  baud.  Some,  or  many,  of  the  tones 
may  not  be  used  (or  equivalently  have  amplitudes  of  zero) 
during  any  or  all  bauds  of  the  packet.  For  example,  to 
generate  bandpass  signals  between  frequencies  fx   and  f2, 
only  tones  between  harmonics  k1=f1/Af  and  k2=f2/Af  will  be 
allowed  non-zero  amplitudes.  Here,  the  maximum  number  of 
tones  will  be  K=k2-k1+l  and  the  signal  bandwidth  will  be 
W=K*Af.  Note  that  the  time  bandwidth  product  of  the  entire 
signal  packet,  TW=LAT*KAf=LK,  is  equal  to  LK,  the  total 
number  of  symbols  that  can  be  sent  in  the  packet. 

Let  us  now  provide  a  mathematical  description  of  the 
signal  packet.  Let  the  1th  baud  of  the  transmit  signal  be 
described  by; 


*i(u)  =  I  xu(u)  (1) 

k 

where, 

xi*(u)    =  Aikcos(27rkAfu   +   0uJ     ;    0<u<AT.  (2) 


Here,  u  is  time  referenced  to  the  beginning  of  the  baud. 
Actual  real  time  is  t=t0+lAT+u  where  t0  is  the  time  of 
initiation  of  the  0th  baud,  i.e.,  the  beginning  time  of  the 
packet.  Now  the  discrete  time  signal  corresponding  to  the 
1th  baud  is  found  by  sampling  (1)  and  (2)  at  the  sampling 
intervals  At=l/fx  and  is  given  by; 

*i(n)  =  I  xlk(n)  (3) 

where, 

*ik(n)  =  AIkcos(27rkn/kx  +  <p^)     ;  0<n<kx-l.   (4) 

Here,  n  is  discrete  time  referenced  to  the  beginning  of 
the  baud.  Note  that,  in  general,  k  may  take  on  all  integer 
values  between  0  and  kx/2.  We  will  refer  to  k  as  the 
"harmonic  number"  of  the  MFM  tone  of  frequency  kAf.  It  is 
useful  to  note  that  a  baud  interval,  i.e.,  a  time  AT, 
contains  exactly  k  cycles  of  tone  k.  Thus,  adjacent  tones 
differ  by  one  in  the  number  of  cycles  they  make  during  a 
baud. 

Consider  the  kx  point  Discrete  Fourier  Transform  (DFT)  of 
(3).  It  is  given  by; 


Xi(k')  =  I  X^Jk')  (5) 

k 


where, 


X^(k')  =  J5kxAIk{exp(-j0li)6(k'-k) 

+  exp(j0li)6(k'-(kx-k))  }  ;  0<k'<kx-l.   (6) 

Thus,  it  is  apparent  that  the  discrete  time  signal  (3)  is 
given  by  the  3^  point  Inverse  Discrete  Fourier  Transform 
(IDFT)  of  (5)  ,  namely; 

xx(n)  =  IDFTtX^k)  ]  (7) 

where  from  (6)  it  is  clear  that, 


XL(k)  =  J^A^expt-jc^)  ;    0<k<kx/2 


X^k.-k)  =  ^A^exptj^)   ;   0<k<kx/2 


(8) 


To  summarize,  the  1th  baud  is  generated  by  taking  the  IDFT 
of  a  complex  valued  array  of  length  kx.  The  first  half  of 
the  array  is  loaded  with  the  amplitudes  and  phases  of  the 
tones  to  be  included  in  the  MFM  signal  at  the  corresponding 
harmonic  numbers.  The  second  half  of  the  array  is  loaded 
with  the  complex  conjugate  of  the  values  in  the  first  half 
of  the  array  at  the  image  harmonics  (the  image  harmonic  of  k 
is  kx-k)  .  This  data  structure  is  shown  for  kx=16  in  Fig.  2. 
The  complex  conjugate  symmetry  about  the  midpoint  insures 
that  the  IDFT,  also  an  array  of  kx  points,  will  be  real 
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valued.  The  values  in  this 
array,  xx(n),  are  the 
discrete  time  samples  of  the 
analog  transmit  signal. 
Clocking  these  samples  out 
thru  a  D/A  converter  (they 
are  stored  as  binary  numbers 
in  the  computer)  at  fx 
samples  per  second  completes 
the  generation  of  the  1th 
baud.  The  total  signal 
packet  is  generated  by  an  L- 
fold  repetition  of  this 
process . 
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Figure  2  Data  Structure  for  IDFT 


Figure  3  Block  Diagram  of  Transmitter, 


A  block  diagram  of  the  transmitter  is  shown  in  Fig.  3. 
At  NPS  we  have  implemented  this  using  an  IBM  PC/AT  as  a  host 
computer.  Samples  may  be  clocked  out  at  any  clock  rate  up  to 
250  KHz;  this  limit  is  set  by  the  maximum  byte  transfer  rate 
of  the  AT"  s  DMA  channel  [ref  5].  The  maximum  packet  size  is 
61440  samples;  this  limit  is  set,  currently,  at  one  sector 
of  the  AT's  memory.  We  have  been  using  this  system  to 
generate  bandpass  signals  in  the  10-14  KHz  band  and  16-20 
KHz  band  using  a  clock  rate  of  61440  samples  per  second.  The 
relevant  signal  parameters  are  given  in  Appendix  1. 

Demodulation  of  MFM  is  accomplished  by  a  process  inverse 
to  its  generation.  Given  the  time  domain  signal  at  the 
receiver  xL(u)  on  the  interval  0<u<AT:  The  signal  is 
sampled  at  fx  samples  per  second  and  converted  to  digital 
format  with  an  A/D  converter.  The  kx  real  values  thus 
obtained  are  loaded  into  a  kx  point  complex  valued  array 
(the  imaginary  parts  are  set  to  zero) .  This  becomes  the 
array  xL(n)  for  the  1th  baud.  Its  K^   point  DFT  yields  the 
complex  valued  array  XL(k)  containing,  in  its  first  half, 
the  amplitude  and  phase  modulation  information,  A^  and  0^, 
of  the  K  harmonics   employed  in  the  generation  of  the 
transmit  signal. 

Notice  that  the  demodulation  operation  is  linear 
(ignoring  for  the  moment  the  non-linear  nature  of  the  A/D 
conversion) .  Also  notice  that  only  K  out  of  kx  complex 
values  are  retained  for  decoding.  The  upper  half  of  the 
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array  is  redundant;  it  contains  the  complex  conjugates  of 
the  values  at  the  image  freguencies.  Also,  non-used 
harmonics  in  the  lower  half  of  the  array  will,  in  the 
absence  of  noise,  contain  zero  values  and  can  be  discarded. 
When  they  do  contain  noise  values,  their  discard  is 
eguivalent  to  filtering  out  the  unused  portion  of  the 
spectrum  at  the  ouput  of  the  demodulation  stage. 

Later  in  this  report,  it  will  be  shown  that  DFT 
demodulation  of  the  signal  is  eguivalent  to   correlation,  or 
matched  filter,  reception  for  each  of  the  K  tones.  It  will 
also  be  shown  that  the  tones  are  all  orthogonal;  in  discrete 
time  on  the  interval  0<n<kx-l,  and  in  continuous  time  on  the 
interval  0<u<AT.  Of  course,  signals  in  different  bauds  are 
orthogonal,  since  they  do  not  overlap  in  (real)  time.  Thus, 
the  K  tones  of  the  L  bauds  form  a  set  of  LK  orthogonal 
signals  so  the  response  of  a  matched  filter,  or  correlator, 
to  any  of  the  tones  (regardless  of  its  modulation)  other 
than  the  one  to  which  it  is  matched  is  zero.  The  fact  that 
the  receiver  is  linear,  and  a  matched  filter  assures  that  it 
is  optimal  for  demodulation  of  MFM,  that  is  it  maximizes 
signal-to-noise  ratio,  in  the  presence  of  additive  white 
noise.  This  is  so  in  spite  of  unknown  link  attenuation  at 
each  of  the  freguencies.  We  note  that  the  channel  need  only 
have  constant  attenuation  over  narrow  bands  corresponding  to 
the  bandwidths  of  the  individual  tones,  Af,  and  not  over  the 
entire  bandwidth  W.  Thus  the  demodulation  is  still  optimal 
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for  a  wideband  signal  propagating  thru  a  channel  with 
substantial  variation  in  gain  across  the  band.  A  block 
diagram  of  the  MFM  receiver  is  shown  in  Fig. 4. 


CUr.t 


IjiUA) 


WD 


±^ 


DFT 


W 


S  i 4  vi  a  \    D €  wi  od  w  \ot'0>v 


T3au,d  i" 


Figure  4  Block  diagram  of  receiver 

Properties  of  MFM  Signals 

In  this  section,  a  number  of  important  properties  of  MFM 
signals  will  be  presented.  In  almost  every  case  we  shall  be 
referring  to  a  single  baud,  say  the  1th  baud,  so,  unless 
required  for  clarity  of  presentation,  the  subscript  1  will 
be  dropped. 


a.  Orthogonality  of  subsignals. 
In  continuous  time; 

'  0   ;  k*i 


AT 


J   xk(u)Xi(u)du  =  < 


^(Ak)2AT   ;  k=i 


(9) 
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and    in   discrete   time; 

!0      ;    k*  i 
(10) 
•2  ( Ak )    kx      ;    k=  1 . 


Eqns  (9)  and  (10)  are  easily  established  using  (2)  and  (4). 

b.  Autocorrelation  function  (acf)  and  power  spectral  density 
(psd)  of  MFM  signals. 

The  circular  (or  periodic)  acf  of  an  MFM  signal  is 
defined  as; 


rx(p)  =  I  x(n)x(n©p)   ;  0<p<kx-l  (11) 


where  ©  signifies  a  left  circular  shift.  Notice  that  a 
circular  shift  of  any  of  the  harmonic  components  of  x(n)  is 
still  a  harmonic  component  of  the  same  frequency  but  with  a 
different  phase  shift.  Combine  this  observation  with  the 
orthogonality  property  (10)  and  (11)  is  easily  evaluated  to 
be  ; 


rx(P)  =  I  rk(p)  (12) 

k 

with, 

rk(p)  =  ^(Ak)2kxcos(27Tkp/kx).  (13) 
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The  acf  is,  as  expected,  independent  of  any  phase  modulation 
on  the  tones;  however,  it  is  modified  by  amplitude 
moduation. 

It  is  well  known  that  [ref  6]    ; 

DFT[rx(p)]  =  X*(k)X(k)  (14) 

and  the  discrete  psd  is  defined  as, 

Sx(k)  =  X*(k)X(k)/kx  (15) 

From  (8)  it  is  seen  that, 

Sx(k)  =  Sx(kx-k)  =  ^,(Ak)2kx  ;   0<k<kx/2.        (16) 

Two  interesting,  and  important  examples,  are  provided  by 
white  lowpass  and  white  bandpass  sequences.  The  white 
lowpass  sequence  is  defined  by; 

I  A   ;  0<k<K-l 
(17) 
0   ;  K<k<kx/2 

Its  acf  is  determined  by  evaluting  (13)  to  be; 


rx(p)    =      ^A2kxcos{7r(K-l)p/kx}    sin(irKp/)<Ly)  (18) 

sin(7rp/kx) 
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The  white  bandpass  sequence  is  defined  by; 

(A      ;  ki<k<k2 
Ak  =  \  (19) 

VO   ;  other  k  between  0  and  kx/2. 

Recalling  that  the  number  of  harmonics  K=k2-k1+l  and 
defining  a  midband  harmonic  number  k0=(k2+k1)/2  then  (13)  can 
also  be  evaluated  in  closed  form.  For  an  even  number  of 
harmonics  in  the  sum  (  K  an  even  number) ,  the  acf  of  the 
white  bandpass  sequence  is  given  by; 


rx(p)  =  ^A2kxcos{27rk0p/kx}  sin(7rKp/kxl         (20) 

sin(?rp/kx) 


Figure  5  shows  the  acf  and  psd  of  a  white  bandpass 
sequence  with  a  kx  of  128,  a  kx  of  47  and  a  k2  of  54  (K=8, 
k0=50.5).  Note  the  nulls  in  the  envelope  of  the  acf  occur  at 
multiples  of  p=kx/K.  In  Fig.  5  this  is  p=16  (as  it  is  for 
all  the  signals  in  Appendix  1) . 

To  summarize,  the  circular  acf  of  MFM  signals  is  given  by 
(13)  for  general  amplitude  modulated  signals,  but  simplifies 
to  (18)  for  white  lowpass  sequences  and  to  (20)  for  white 
bandpass  sequences. 

For  test  purposes,  or  synchronization,  we  may  elect  to 
repeat  a  signal  for  several  bauds.  Suppose  x(n)  is  repeated 
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Figure  5  ACF  and  PSD  of  a  white  bandpass  sequence 
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L1  times;   it  will  have  L'k,.  samples  points.  Its  circular 
acf  is  given  by; 


rx(p)  =  I  x(n)x(n©p)  ;  0<p<h'K-^-  (21) 


but  the  harmonics  are  still  orthogonal  over  the  repeated 
interval  so  that; 


rx(p)  =  I  rk(p)   ;   0<p<L'kx-l  (22) 

k 

with, 

rk(p)  =  ^(Ak)2L'kxcos(27rkp/kx)  .  (23) 


Note  that  the  acf  has  L'  times  the  amplitude  of  the  acf  of  a 
single  baud,  and  since  it  is  still  periodic  with  period  kx, 
there  are  peaks  every  kx  points.  For  signals  with  K  constant 
amplitude  harmonics,  such  as  the  white  sequences  described 
above,  there  will  be  periodic  peaks  of  amplitude; 

rx(mkx)  =  A2L'kxK/2   ;   0<m<L'-l  (24) 

An  acf  and  psd  for  the  signal  of  Fig.  5  repeated  four  times 
is  shown  in  Fig.  6. 

c.  Matched  filter  (correlator)  detection  of  MFM. 

MFM  signals  may  be  further  resolved  into  their  quadrature 
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Figure  6  Acf  and  psd  of  baud  repeated  four  times 


18 


components.  That  is; 

xk(n)  =  xIk(n)  +  xQk(n)  (25) 

where, 

xIk(n)  =  Re[X(k)]cos(27rkn/kx)/kx 
and,  (26) 

xQk(n)  =  -lm[X(k)  ]sin(27rkn/kx)/kx 

with    X(k)  =  -2kxAkexp( j0k)  the  complex  DFT  coefficient  of 
the  kth  harmonic.  Note  that  xIk(n)  and  xQk(n)  are  orthogonal 
over  a  baud  interval  (or  intervals)  and  they  are  orthogonal 
to  the  in-phase  and  quadrature  components  of  all  other 
harmonics  too.  So,  when  resolved  into  quadrature  components, 
the  MFM  signal  space  actually  contains  2K  orthogonal  signals 
in  each  baud  or  2KL  orthogonal  signals  in  a  packet. 
The  impulse  response  for  a  filter  matched  to  the  ktb  tone  of 
the  1th  baud  (we  are  still  supressing  1  in  our  notation) 
consists  in  fact  of  a  pair  of  filters;  one  filter  is  matched 
to  the  in-phase  signal  and  one  to  the  quadrature  signal. 
Their  impulse  responses  are  given  by; 

hIk(n)  =  ClXlk(kx-n)   ;    0<n<kx-l  (27) 

for  the  in-phase  component  and, 

hQk(n)  =  CQxQk(kx-n)   ;    0<n<kx-l  (28) 
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for  the  quadrature  component.  CT   and  CQ  are  arbitrary 
constants.  The  output  of  the  discrete  time  matched  filter 
pair,  when  excited  with  x(n),  the  1th  signal  baud,  is  given 
by  the  convolutions; 


Yik(n)  =  hIk(n)*x(n) 
and,  (29) 

yQk(n)  =  hQk(n)*x(n)  . 


Evaluating  these  at  the  end  of  the  baud,  i.e.,  at  n=kx,  and 
choosing  Cz   and  CQ  to  make  the  filters  independent  of  the 
modulation  yields; 


Yik(kx)  =  I  x(n)cos(27rkn/kx)=  Re[X(k)] 
and,  (3  0) 

yQk(kx)  =  I   x(n)sin(2?rkn/kx)  =  -Im[X(k)]. 
n 


Thus,  the  real  and  imaginary  components  of  the  DFT  of 
x(n)  at  frequency  k  are  identical  to  the  outputs  of  an  in- 
phase  and  quadrature  pair  of  filters,  matched  to  frequency  k 
and  sampled  at  the  end  of  the  baud.  Put  another  way,  the  DFT 
of  x(n)  gives  us  K  complex  outputs,  or  2K  real  outputs, 
which  are  equivalent  to  the  outputs  of  K  pairs  of  matched 
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filters,  each  matched  to  one  of  the  signals  from  our 
orthogonal  set  of  2K  signals,  sampled  at  the  end  of  the 
signal.  The  linearity  of  the  DFT  insures  that  the  response 
to  signal  plus  noise  of  the  DFT  is  the  same  as  the  response 
of  the  matched  filter.  In  particular,  the  DFT  maximizes  the 
signal-to-noise  ratio  for  detecting  the  amplitudes  of  the 
guadrature  components  of  all  the  harmonics  in  the  MFM 
signal.  In  this  sense,  the  DFT  is  the  optimum  reciever  for 
MFM. 

III.  PERFORMANCE  OF  MFM 
Signal  Received  Plus  Noise 

In  general,  the  signal  arrives  at  the  receiver  over  a 
link  that  will  introduce  freguency  dependent  attenuation  and 
phase  shifts.  It  may  introduce  a  significant  bulk  delay 
(  for  example  a  fraction  of  a  second  in  satellite 
communications   and  perhaps  several  seconds  in  acoustic 
communications) .  The  bulk  delay  may  vary  with  time  due  to 
transmitter  and/or  receiver  dynamics  introducing  time  and 
freguency  dilation  (generalized  doppler  shift) .  The  transfer 
characteristics  of  the  channel  may  be  random  in  some  way  as 
in  channels  with  fading.  Inevitably  noise  will  be  added  to 
the  signal,  either  from  the  environment,  the  electronics  or 
both. 

In  this  report  it  will  be  assumed  that  the  channel 
transfer  function  consists  only  of  known  attenuation  and 
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that  the  noise  is  additive  white  guassian  noise.  In  such  a 
channel,  synchronization  may  be  assumed  and  coherent 
reception  of  the  in-phase  and  quadrature  symbols  is  possible 
as  described  in  the  previous  section.  The  very  important 
problems  that  delay  and  doppler  shift  create  for 
synchronization  and  the  performance  of  synch  algorithims  are 
being  treated  in  a  separate  report.  (Multipath  and  fading 
problems  in  MFM  are  the  subject  of  a  future  study.  They  are 
not  considered  to  be  serious  problems  in  the  near  vertical 
acoustic  link  in  whose  development  we  are  currently 
involved.  However,  synchronization  is  a  principal  researh 
issue  because  of  the  lengthy  acoustic  propagation  times.) 
Since  the  signals  are  being  processed  baud-wise,  we  shall 
continue  to  suppress  the  baud  notation,  understanding  that 
our  results  apply  to  the  "1th  baud"  for  arbitrary  1.  Baud 
identity  will  be  reintroduced  only  as  required  for  clarity. 

With  synchronization,  the  recieved  signal  plus  noise 
will  be  sampled  at  fx  samples  per  second,  and  the  first 
sample  will  be  triggered  at  time  u=0.  This  sampling  will 
yield  the  discrete  time  signal; 


y(n)  =  X  (2Pk)%cos(27rkn/kx  +  0k)   +   w(n)  ;  0<n<kx-l   (31) 
k 


where  Pk  is  the  received  signal  power  of  the  kth  tone  and 
w(n)  is  a  sequence  of  white  noise  samples.   Let  us  relate 
the  statistical  properties  of  the  white  noise  sequence  w(n) 
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to  those  of  the  analog  continuous  time  white  noise  signal 
w(u)  . 

Assume  that  the  received  data  y(u)  has  been  ideally 
bandlimited  to   one  half  the  sampling  frequency  such  that 
there  is  no  power  in  frequencies  greater  than  or  equal  to 
fx/2.  If  w(u)  has  a  white  PSD  equal  to  N0/2  then; 

E[w(u) ]  =  E[w(n) ]  =  0  (32) 

and, 

a2    =   Var[w(u)]  =  N0fx/2  =  Var[w(n)].  (33) 

Also,  sampling  at  twice  the  highest  frequency  assures  that 
the  noise  samples  of  white  noise  are  uncorrelated,  i.e., 

E[w(n)w(m)]  =0   ;  n  h  m.  (34) 


Signal-to-Noise  Ratios  (SNR's)  and  the  DFT 

Let  the  receiver  input  SNR  be  defined  as  the  signal  power 
in  bandwidth  W  divided  by  the  noise  power  in  bandwidth  W. 
Also  define  the  average  tone  signal  power  as; 


Po  =  I    Pk/K  (35) 

K 

then, 

SNR!    =    KP0/WN0    =    P0/AfN0    =    P0kx/(2a2)     .  (36) 
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Note  that  the  narrowband  input  SNR; 

SNRNB  =  Pk/AfN0  =  Pkkx/(2a2)  (37) 

is  the  same  as  the  wideband  SNR  of  (36)  if  the  tone  power  is 
the  same  as  the  average  power,  a  certainty  if  all  tones  have 
equal  received  powers  (Not  necessarily  the  same  as  all  tones 
having  the  same  transmit  powers  if  different  tones  suffer 
different  attenuations.  If  the  channel  frequency  response  is 
known  at  least  approximately,  then  it  may  be  feasible  to 
pre-emphasize  the  transmit  tone  powers  such  that  the 
received  signal  is  approximately  white.  Of  course,  this  may 
not  be  practical  if  there  are  large  nulls  in  the  frequency 
response. ) 

Now  consider  the  k^-point  DFT  coefficients  of  the  input 
sequence; 

Y(k)  =   S(k)  +  W(k)  (38) 

where, 

S(k)  =  ^(2Pk)%kxexp(J0k)  ;  kx<k<k2  (39) 

and  the  W(k)  are  the  DFT  coefficients  of  the  white  noise 
sequence.  It  is  easily  shown  that  the  W(k)  have  the 
following  statistical  properties; 

E[W(k) ]  =  0  (40) 
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E[Re{W(k)  }2  ]    =    E[Im{W(k)}2]    =    k^a2  /2  (41) 

E[Re{W(k) }Im{W(k) }]    =    0.  (42) 

Furthermore,  the  Real  and  Imaginary  parts  of  W(k)  and  W(i) 
are  all  uncorrelated  with  one  another  for  k  ^  i. 

To  summarize,  the  2K  coefficients  given  by  the  Real  and 
Imaginary  parts  of  the  K  points  of  the  DFT  corresponding  to 
the  tones  used  in  the  MFM  transmit  signal  have  means  given 
by  the  Real  and  Imaginary  parts  of  (39)  and  variances  given 
by  (40).  They  are  uncorrelated  random  variables.  If  the 
input  noise  seguence  w(n)  is  gaussian,  they  are  uncorrelated 
gaussian  random  variables. 

Define  the  output  SNR  of  receiver  as  the  ratio  of  the 
sguare  of  the  mean  to  the  variance  for  each  of  these  2K 
coefficients.  These  ratios  are  the  maximum  since,  as  we  have 
shown  previously,  the  DFT  coefficients  are  identical  to  the 
output  of  matched  filters.  Thus, 

[SNRI]k  =  Pkkx(cos0k)2/a2 
and,  (43) 

[SNRQ]k  =  Pkkx(sin0k)2/a2 

It  is  useful  to  note  from  (37)  that; 
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P^K/O2     =    2SNRNB.  (44) 

For  the  purpose  of  analysis,  define  the  unit  variance, 
uncorrelated  random  variables; 


Ik  =  Re[Y(k)]/{Var(Re[Y(k)])}% 
and,  (45) 

Qk  =  Im[Y(k)]/{Var(Im[Y(k)]) }% 


and  note  from  (43)  that; 

E[Ik]  =  ±{[SNRI]k}%  =  (P,kya!  }%cos0k 
and,  (46) 

E[Qk]  =  ±{[SNRQ]k}%  =  {Pkkx/a2  }%sin0k 

where  the  ±  depends  on  the  phase  0k.  Furthermore,  if  the 
noise  is  gaussian,  these  random  variables  are  unit  variance, 
statistically  independent  gaussian  random  variables  with 
mean  values  given  by  (46) . 

Bit  Error  Probabilities 

In  this  section  results  are  presented  for  bit  error 
probabilities  of  three  types  of  MFM:  MFBPSK,  MFQPSK  and 
MF16-QAM.  These  modulation  formats  carry  one,  two  and  four 
bits  per  tone  per  baud.  Thus  they  carry  one,  two  and  four 
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Figure  7  MFM  Phase  Modualation  Formats  (Hex) 

bits  per  Hz  of  bandwidth  per  second  or  equivalently  KL,  2KL 
and  4KL  bits  per  packet  respectively.  The  Ik  and  Qk 
variables,  in  the  absence  of  noise,  are  shown  in  Fig.  7   for 
each  of  the  three  modulation  formats.  The  magnitudes  of  the 
tones  have  been  choosen  for  equal  tone  powers  in  MFBPSK  and 
MFQPSK  and  for  an  equal  average  tone  power  for  MF16-QAM 
assuming  each  of  the  16  symbols  is  equally  probable.   This 
makes  it  possible  to  compare  modulation  types  under  the 
condition  of  equal  total  (or  average  total)  transmit  power. 


a.  MFBPSK 

In  MFBPSK,  the  quadrature  output  is  not  used.  Assuming 
0=0,  then  a  correct  bit  detection  is  made  anytime  I  >  0, 
that  is,  anywhere  in  the  right  half  plane  of  Fig.  7a.  The 
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probability  of  bit  error  is  given  by; 


PE  =  1  -  Pr[I  >  0|  0=0]  =  Q{E[I]}  (47) 


where  from  (44)  and  (46), 


E[I]  =  {2SNRNB}*  (48) 


and  Q(x)  is  the  Q-function  [ref  1].  The  function  is  bounded 
above  by  the  exponential; 


Q(x)  <  [exp(-xV2)]/(27r)%  (49) 


and  for  x  greater  than  two,  the  bound  is  a  very  good 
approximation  to  Q(x). 

Eqn  (47)  is  plotted  in  Fig.  8  as  a  function  of  SNRNB. 

b.  MFQPSK 

In  MFQPSK,  both  the  in  phase  and  quadrature  channels 
carry  a  bit  of  information.  However,  since  they  are 
statistically  independent,  we  can  decode  them  independently. 
For  example,  if  0  =  tt/4   corresponds  to  logical  00  and 

0  =  -7r/4  corresponds  to  logical  10,  then  the  least 
significant  bit  will  be  correctly  decoded  as  a  0  so  long  as 
I  >  0  when  either  00  or  10  is  sent.  Thus,  the  probability 
of  bit  error  is  given  by; 
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PE   =    1    -    Pr[I    >    0|    <p   =   ±7T/4]    =   Q{E[I]}  (50) 


where   again    from    (44)    and    (46) , 


E[I]    =    (SNRNB}%  (51) 


Eqn    (49)    is   also   shown    in   Fig   8. 

c.    MF16-QAM 

In  MF16-QAM,  the  in  phase  and  quadrature  channels  both 
carry  two  bits  of  information.  Again,  since  they  are 
statistically  independent,  we  can  decode  them  independently 
and  by  symmetry,  the  bit  error  probabilities  calculated  for 
one  channel  will  be  true  for  the  other  channel  as  well. 

Let  us  concentrate  on  the  in  phase  channel.  There  are 
four  symbols;  two  with  amplitudes  ±(2PL)%  and  two  with 
amplitudes  ±(2PH)%.  We  want  the  average  power  in  the  in 
phase  channel  to  be  Pk/2,  the  same  as  QPSK.  Assuming  that 
low  and  high  power  symbols  are  equally  probable; 

hPL   +  hPE   =   Pk/2  (52) 

and    for    16-QAM  we  make    (2PH)%  =   3(2?^*   then, 

PL   =   Pk/10   and   PH  =   9Pk/10.    Thus,    referring   to    (46)    we   see 

that    ; 

29 


E[I|PL]    =   E[Re{Y(k)  |PL]/{Var[Re{Y(k)}]}%   =    (PLkx/a2)%       (53) 
or, 

E[I|PL]    =    ±{SNRNB/5}%      =    ±3  (54) 

Similarly; 

E[I|PH]  =  ±{9SNRNB/5}*  =  ±3a  (55) 


Now  referring  to  Fig.  7,  assuming  equally  likely  symbols, 
the  probability  of  correct  decoding  of  in  phase   symbols  is 
given  by; 

Pc  =  h(    Pr{0<I<2a|PL}  +  Pr{2a<l|PH)  )  =  1  -  3Q(a)/2    (56) 

Since  each  symbol  contains  2  bits,  then  we  can  define  an 
equivalent  bit  error  probability  such  that  two  bit  sequences 
will  have  the  same  probability  of  correct  decoding.  That  is 
the  probability  of  correctly  decoding  two  bit  sequences  of 
indepedent  bits  is  given  by; 

Pc  =  (1  "  PEeq)2  •  (5V) 

Equating  (55)  and  (56)  gives  an  equivalent  bit  error 
probability  for  MF16-QAM; 
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where, 


PEeq   =    1    -    {1    -    3Q(a)/2}* 


a    =     |E[l|PJ  |     =    {SNRNB/5}*. 


(58) 


(59) 


The  equivalent  probability  of  bit  error  for  MF16-QAM,  (58) 
is  shown  in  Fig.  8  as  a  function  of  SNRNB. 


Figure  8  MFM  Bit  Error  Probabilites  vs.  SNR 


tJB 
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Optimization 

Bit  error  probabilities  are  directly  reduced  by 
increasing  SNRNB  for  MFM  as  illustrated  clearly  in  the  three 
types  of  modulation  described  above.  From  (37) ; 

SNRNB  =  Pk/AfN0  =  PkAT/N0  (60) 

with, 

AT  =  T/L.  (61) 

Now  given  a  fixed  power  for  the  tones  and  fixed  additive 
noise  level,  then  we  want  to  maximize  AT  or  minimize  L  and 
Af.  What  restricts  us  from  choosing  L=l  and  AT=T?  Only  the 
coherence  time  of  the  channel.  If  the  channel  is  stationary 
over  times  greater  than  the  packet  length  T,  then  the  packet 
should  consist  of  a  single  baud  (L=l)  of  length  T,  tones 
will  be  spaced  at  1/T  Hz  and  there  will  be  K=TW  tones  in  the 
channel  bandwidth  W.  If  the  coherence  time  of  the  channel  is 
less  than  the  packet  length  then  AT<TC  and  there  must  be 
L>T/TC  bauds  in  the  packet.  As  an  example,  for  the  near 
vertical  acoustic  link,  it  is  expected  that  Tc  will  be  about 
.025  seconds  and  packets  are  of  length  T=.066  seconds.  Thus, 
packets  must  be  made  up  of  at  least  3  bauds  and  tone  spacing 
Af  must  be  no  less  than  40  Hz. 

As  another  example,  suppose  we  are  using  10  bauds  in  a 
packet  and  achieving  a  SNRNB  of  ten  allowing  transmission  of 
data  at,  10"4'7,  a  satifactorily  low  BER  using  MFBPSK.  This 
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modulation  provides  a  data  rate  of  one  bit  per  Hz  of  channel 
bandwidth.  Now  suppose  the  channel  is  sufficiently  stable 
that  it  is  possible  to  use  only  one  baud  per  packet,  but  ten 
times  as  many  tones.  We  can  now  expect  an  SNRNB  of  100  and 
obtain  the  same  BER  using  MF16-QAM.  This  modulation  gives  a 
data  rate  of  four  bits  per  Hz  of  channel  bandwidth,  that  is 
four  times  as  great  as  the  MFM  using  shorter  bauds. 

Redundancy 

Redundancy  may  be  introduced  in  frequency,  time,  or  both 
in  order  to  improve  SNRNB  and  reduce  bit  error  probability. 
This  is  done  at  the  expense  of  bit  rate.  Recall  that  there 
are  2KL  orthogonal  signals  in  a  packet.  If  an  information 
symbol  is  repeated  M  times  on  M  of  the  2KL  signals  choosen 
in  a  known  pattern,  the  pattern  may  be  psuedo-random  for 
example,  then  the  M  outputs  of  the  2KL  DFT  coefficients  can 
be  added  prior  to  decoding.  This  assumes  that  coherence  is 
maintained  amoung  the  M  signals  in  transmission.  Given 
coherence,  then  SNRNB  will  be  increased  by  M^.  Of  course  the 
bit  rate  will  be  reduced  by  a  factor  of  M.  This  is  similar 
to  the  effect  introduced  by  1/M  rate  convolutional  codes 
[ref  7] . 

More  will  be  presented  on  coding  and  on  signal  spreading 
techniques  using  MFM  in  a  subsequent  report.  It  is  worthy  of 
note  at  this  point  that  redundancy  can  be  useful  to  combat 
particular  channel  problems  like  fading  or  burst  noise. 
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IV.  DISCUSSION  AND  CONCLUSIONS 
MFM  is  a  new  digital  signal  processing  oriented 
communications  signal  ideally  suited  for  computer  to 
computer  links  and  computer  based  information  exchange 
networks.  In  this  report  we  have  presented  the  theory  behind 
MFM  signal  generation  and  reception,  that  is  modulation  and 
demodulation,  using  the  host  computer's  DSP  algorithims  and 
hardware. 

Spectral  and  temporal  properties  of  MFM  have  been  derived 
and  examples  presented  for  bandpass  and  lowpass  channels. 
Finally  MFM  performance,  in  terms  of  bit  error  rates,  has 
been  calculated  for  MFBPSK,  MFQPSK  and  MF16-QAM.  These 
signals  provide,  respectively,  one,  two  and  four  bits  per 
second  per  Hz  of  channel  bandwidth.  We  have  also  seen  that 
for  channels  with  attenuation  only,  optimum  performance, 
that  is  minimum  error  rate  for  a  given  data  rate  or  maximum 
data  rate  for  a  given  error  rate,  is  obtained  by  using  the 
minimum  tone  spacing  and  hence  the  maximum  number  of  tones 
in  the  available  bandwidth.  The  minimum  tone  spacing  is  the 
reciprocal  of  the  packet  length,  or  the  coherence  time  of 
the  channel,  whichever  is  smallest. 

In  many  channels,  particuarily  those  where  propagation 
over  an  unknown  and/or  time  varying  path  is  involved,  there 
is  uncertainty  in  the  time  delay  and  doppler  dilation 
introduced  by  the  propagation.  This  introduces  the  problem 
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of  synchronization.  In  order  to  have  coherent  reception  of 
MFM,  synchronization  must  be  obtained  with  respect  to  time 
of  arrival  of  each  baud  and  with  respect  to  sampling 
frequency  for  the  analog  waveform,  a  waveform  that  may  have 
been  dilated  by  the  channel.  Results  on  synchronization 
error  and  sampling  frequency  error  as  they  effect  the 
performance  of  MFM  will  be  presented  in  a  separate  NPS 
Technical  Report.  Synchronization  algorithims,  as  well  as 
differential  MFM,  a  less  synch  sensitive  form  of  MFM  will  be 
discussed  there  too. 


35 


REFERENCES 


1.  Couch,  L.  W.  II,  Digital  and  Analog  Communications 
Systems,  2nd  ed.,  Macmillan,  1987. 

2.  Ricci,  F.  J.  &  Schutzer,  D.  ,  U.S.  Military 
Communications,  Computer  Science  Press,  1986. 

3.  Kalet,  Irving,  "The  Multitone  Channel",  IEEE  Trans,  on 
Communications,  vol  37,  no.  2,  Feb  1989,  pp  119-124. 

4.  Moose,  P.  H. ,  "Submarine  Acoustic  Tactical  Data  Link", 
Proc.  of  MILCOM  86.,  Oct  1986,  Monterey,  CA. 

5.  Childs,  Robert  Daniel,  "High  Speed  Output  Interface  for  a 
Multifrequency  Quatenary  Phase  Shift  Keyed  Signal  Generated 
on  an  Industry  Standard  Computer",  MSEE  Thesis,  Dec  1988, 
Naval  Postgraduate  School,  Monterey,  CA. 

6.  Strum,  Robert  D. ,  and  Kirk,  Donald  E.,  First  Principles 
of  Discrete  Signal  Processing,  Addison-Wesley,  1988. 

7.  Lin,  S.,  and  Costello,  D.  J.  Jr.,  Error  Control  Coding, 
Prentice  Hall,  1983. 


36 


APPENDIX  1 
DESIGN  PARAMETERS  FOR  1/1 5TH  SECOND 
SIGNAL  PACKET  IN  A 
16-2 OKHZ  BANDPASS  CHANNEL 


AT  sec 

1/240 

1/120 

1/60 

1/30 

1/15 

L 

16 

8 

4 

2 

1 

Af 

240 

120 

60 

30 

15 

*1 

68 

135 

269 

537 

1073 

*1 

16320 

16200 

16140 

16110 

16095 

k2 

83 

166 

332 

664 

1328 

f2 

19920 

19920 

19920 

19920 

19920 

K 

256 

512 

1024 

2048 

4096 

fv 

61440 

61440 

61440 

61440 

61440 
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